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Abstract 

We examine the effect of the transverse breathing mode on the longitudinal sound propagation 
in a Bose-Einstein condensate in a one-dimensional optical lattice. In particular, we discuss how 
the coupling with the transverse breathing mode influences the sound velocity in an optical lattice. 
Using a variational approach we calculate the dispersion relations for the longitudinal sound mode 
and the transverse breathing mode analytically and find that the shift in the sound velocity from the 
uncoupled result can be large enough to be experimentally relevant. We also find that the effective 
mass of the transverse breathing mode is affected considerably by the coupling to longitudinal 
sound. 
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I. INTRODUCTION 



Repulsive contact interactions between atoms in a dilute Bose gas cause a qualitative 
change in the low-energy excitations of the system. In particular, the noninteracting single- 
particle dispersion equal to h 2 k 2 /2m is replaced with the sound dispersion ck that is linear 
in the momentum The existence of this sound mode is important as it allows for the 
possibility of superfluid flow . Sound propagation in a trapped Bose- Einstein condensate 
was studied experimentally by Andrews et al. jj| 0] and Stamper-Kurn et al. in the 
regime when the condensate can be considered as being homogeneous in one direction. 
The experimental observations were in good agreement with theoretical predictions in that 
case [EQIl. 

Bose-Einstein condensates can also be placed in an optical lattice. The periodic structure 
of the lattice has enabled studies of Bloch oscillations number squeezing [rJ, collapses 
and revivals 11 1, and a number of superfluid properties of Bose gases 1^, 3|- Importantly, 



the Bose-Einstein condensates in an optical lattice realize the Bose-Hubbard model and 
can be used to study the quantum phase transition from the superfluid into a Mott-insulator 
state 0, H El . 

Apart from these phenomena, sound propagation is also possible in a Bose-Einstein con- 
densate in an optical lattice. There exists several theoretical calculations for the sound veloc 
ity of a Bose-Einstein condensate in a one-dimensional optical lattice rJ 12, HI 3 III 21 



All these theories for the sound mode are based on the idea that the condensate wave func- 
tion in site n is, in tight-binding approximation, essentially one dimensional and has the 
form & n (x,y,t) = vp(x, y)\J ' N n {t)e lUn ^ . This implies that the only dynamical variables of 
these theories are the total number of atoms N n (t) and the global phase v n {t) in every site, 
but there is no time dependence in the wave function describing the transverse directions. 
Physically, it is not a priori clear that this is a valid assumption because in the presence of 
repulsive interactions between the atoms the profile ip(x,y) does not remain constant as the 
number of particles in the site is changed. For example, increasing the number of atoms from 
its equilibrium value causes the condensate to expand whereas a reduction of the number of 
atoms causes the condensate to contract. Therefore, number fluctuations in every site of the 
lattice are coupled to collective modes of the condensate in each site. In a one-dimensional 
optical lattice where the lattice runs in the longitudinal direction (the z-direction), the cou- 
pling predominantly excites the breathing mode in the transverse direction (the xy-p\a.ne). 



These transverse modes were studied previously in Ref. [21[, in the approximation that the 
coupling with longitudinal sound can be neglected. While a coupling definitely exists be- 
tween the sound mode and the transverse breathing mode, a coupling between the transverse 
quadrupole modes and the sound mode does not exist because these modes are orthogonal. 
Therefore, the results in Ref. j2l| for the quadrupole modes are not affected by the coupling 
with the sound mode. 

In this paper we go beyond the approximation used in Refs. [3, 17, [r| [r| H^, 21 1 , by 
providing a unified theory of longitudinal sound and the transverse breathing mode in a 
one-dimensional optical lattice. Using a variational approach we obtain the magnitude of 
the change to the sound velocity when the coupling with the transverse breathing mode 
is taken into account. It turns out that the relative shift in the sound velocity from the 
uncoupled result approaches a constant a/3/4 — 1 m —0.13 as the strength of the interactions 
is increased, i.e., in the Thomas-Fermi limit. Therefore, the shift is rather large and should 
be taken into account when quantitative results are required. Likewise the shift in the 
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effective mass of the transverse breathing mode from the result in Ref. j21| also turns out 
to be large and should also be taken into account. 

Intriguingly, our approach leads to a very simple theory of Josephson oscillations coupled 
to just one transverse degree of freedom. Such a theory has some analogies with the theories 



describing damped relative phase dynamics of weakly- coupled condensates |22|,|23[. In these 
theories the damping of the relative phase is caused by the coupling into quasiparticles. Our 
minimal theory can be used to study the phase dynamics of the Bose-Einstein condensate 
in a one dimensional optical lattice and with considerable modifications even the (quasi) 
irreversible damped phase dynamics. 

This paper is organized as follows. In Sec. |Hj we outline the theoretical foundation of our 
theory and explain the variational ansatz we use. In Sec. IIHI we write down the equations 
of motion and solve for the eigenmodes. We then proceed to dicuss, in subsections IIII Al 
and 1111 Bl the dispersion relations for the sound and the transverse breathing mode sepa- 
rately. Previous experimental studies of the condensate sound mode 0, 0] tracked the time 
evolution of a density dip. Therefore in Sec. II VI we apply our theory to solve numerically 
the time evolution of such a dip in the condensate density in an optical lattice. In addition, 
in Sec. II VI we also solve the time evolution of the system that was prepared with a localized 



breathing-mode disturbance. We conclude with a brief discussion of our results in Sec. |V| 



II. VARIATIONAL ANSATZ AND LAGRANGIAN 



The theory we use is similar to one presented elsewhere |21|. In this paper the only dif- 
ference with the theory presented in Ref. jljj is the use of the grand-canonical Hamiltonian, 
i.e., inclusion of the chemical potential term — /j,N into the energy functional. Due to the 
presence of the global phase factors and the atom-number fluctuations in different sites, this 
turns out to be more convenient for our purposes. 

We consider a Bose-Einstein condensate trapped by a harmonic trap with a radial trap- 
ping frequency uj r . The longitudinal trapping frequency u> z <C uj t is assumed to be so small 
as to be irrelevant. The Bose-Einstein condensate also experiences a one-dimensional opti- 
cal lattice in the longitudinal direction and this lattice splits the condensate into a stack of 
weakly-coupled two-dimensional condensates. Furthermore, we use trap units, i.e., the unit 
of energy is huo ri the unit of time is l/ui r , and the unit of length is l r = ^h/muj r , where m 
is the atomic mass. The energy functional for the stack of two-dimensional condensates is 
then 

\{* 2 + y 2 ) + (i) 



+ %|$ n (^)| 2 -/i 



<n,m> 



where the lattice sites are labeled by n and (n, m) indicates nearest neighbours. This energy 
functional is characterized by two tunable paremeters: the strength of the interaction U2D 
in every two-dimensional Bose-Einstein condensate and the strength J of the Josephson 
coupling between the condensates in the neighbouring sites. In terms of the scattering length 
a, the trapping frequency u L in every site due to the optical lattice, and the characteristic 
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length scale II = ^/h/mujL, the interaction strength is given by 



U 2D = MI- Ij- 



(2) 



Moreover, approximating the lattice potential near its maximum by an upside-down 
parabolic potential we obtain for the Josephson coupling 



J 



8tt 2 



I0 r 
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T 
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-(A/4^) 2 



(3) 



where A is the wavelength of the laser beams creating the optical lattice. 

We use a variational approach to study the sound and the transverse breathing modes of 
the condensate. For this purpose we use the gaussian ansatz 



$ n (x,y,t) = 



N(B + e'(t))(l + 5 n (t)) 



exp 



(Bo + e n (t)) (x 2 + 



+ ivn(t)) (4) 



for the two-dimensional wave function at site n. The variational parameters e n = e' n + ie" n 
are complex and describe the amplitude of the breathing mode, whereas 5 n and u n are the 
relative number fluctuation and the global phase at site n, respectively. Furthermore, N 
is the equilibrium number of atoms at the site and Bq gives the equilibrium size of the 
two-dimensional condensate. The latter is obtained by minimizing the equilibrium energy 
functional. 

Using this ansatz we can calculate the energy functional. We find that in equilibrium the 
condensate size parameter B is given by 

Bo = , 1 , (5) 

where, for convenience, we defined the strength of the interaction as 



U 



N 



In addition, the chemical potential is given by 

3 B 



// 



2£?n 



-2J, 



(6) 



(7) 



because this choice for the chemical potential removes the terms linear in 5 n from the energy 
functional. 

As we are interested in the collective modes we must expand the Lagrangian L = T — E 
up to second order in the variational parameters. Due to the technical simplicity of the 
number and global phase fluctuations we choose, for the time being, to treat them exactly 
and only expand in e n . It then turns out that the contribution from the time- derivative term 
in the Lagrangian is 
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whereas the energy functional becomes 



E/N = Hj + J2 



1 + Sr, 
2Bn 



+ UB (l+8 n y{l + ^-\-fiN(l + 8 n 
In Eq. (jOJ) we have explicitly split off the contribution from the Josephson coupling 



Hj = - J ^ l COS ( U n ~ V m) I'nm ~ sin { V n ~ V m ) I, 



" 1 



(9) 



(10) 



<n,m> 



where I mn = J d 2 r§* m (x,y)§ n (x,y) is the overlap integral. Expanding the overlap integral 
up to second order, now also in the number and global phase fluctuations, we readily obtain 
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J 
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(11) 



From this expression it is clear how the global phase fluctuations are coupled to the transverse 
breathing mode. Technically this is due to the imaginary part of the overlap integral which 
has a contribution linear in (e" — e^J. Such a contribution does not exist, for example, for 
the quadrupole modes which preserve the condensate volume. Therefore, quadrupole modes 
will not be affected by the presence of the sound mode. 

In the above result for the Josephson energy Hj we expanded also in terms of the number 
S n and global phase v n fluctuations. In the problem we are focusing on here these fluctuation 
are small and the above procedure is justified. In principle, however, number and global 
phase fluctuations can also be included exactly. By including them exactly we can also 
cap ture the physics of modulational and dynamical instabilities in an optical lattice ji^, HE 
261 ] . but this is outside the scope of the present paper. 



III. EQUATIONS OF MOTION AND THE EIGENMODES 



The results obtained in the previous section enable us to derive the linearized Euler- 
Lagrange equations of motion for the variational parameters. They read 



de' 

(n,m) 



dt 



B> 



-e' n - 2UB 2 5 n — J ^2 (e' n - e' n 



{n,m} 



dt 



J 



{l + 3U)e' n -3B U5 n --J2 



(n,m) 



el - el 



Bn 



(12) 

(13) 
(14) 
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and 

§ = 23 £ [(„„-„ J) (15) 
These equations of motion can also be written in a more formal way as 



dt u \de'J u \du, 



at » V*;J • 

^ ^ BoflfV (18) 



and 

dt dv n 1 j 

While the first form of the equations is eventually needed in the actual calculations, the 
second way of writing the equations provides some additional insight. In particular, by 
inspecting Eqs. ([16| ) -([18 )) it becomes clear how the global phases influence the behaviour of 
the breathing mode through the term Bq {dE/dv n ) in Eq. (j!6|) and how the breathing mode 
influences the dynamics of the global phase through the term —Bo(dE/de' n ) in Eq. (|18j). 
Without these terms we could treat the condensate density fluctuations independently from 
the transverse breathing mode. 

The four first-order differential equations Eqs. ()12 )) -([15 )) for the variational parameters 
can be cast into two coupled second-order differential equations for e' n and 5 n . We are looking 
for solutions of the type e' n = e' k {t) sin (nkX/2) and S n = 5k(t) sin (nkX/2). By inserting these 
into the equations of motion we obtain 



(2£? + 2J(k)) (J- + 2J(A0J e' k + {B 2 Q - l) (2B + 2J(k)) 8 k 



and 



dt 2 

§ - «w ( B « - i - Jik) ) 4 - ^ ( 3 - 3B » - i) 4 (21) 

where J(k) = J (1 — cos(A;A/2)). The two eigensolutions of these equations correspond to 
the longitudinal sound mode and the transverse breathing mode. These equations of motion 
and their solutions are the main results of this paper. 

In Fig. ^ we show a typical behaviour of both the sound mode and the transverse 
breathing-mode frequency as a function of momentum k for two different values of the in- 
teraction strength. This figure demonstrates how the small k behaviour of the sound modes 
is linear whereas that of the breathing mode is quadratic. It also shows that, both modes 
become stiffer as the strength of the interactions is increased. We devote the following two 
subsections to the discussion of both modes separately. 



(20) 
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FIG. 1: Eigenmode frequencies of the Eqs. (|20|) and (|21jl as a function of momentum, when J = 0.1, 
U = 100 (solid line), and U = 10 (dashed line). 



A. The longitudinal sound mode 

Theories of bosons in a lattice that reduce the problem to that of a global phase and an 
atom number in each site, can be used to solve the dispersion relation u>s{k) of the sound 
mode. We refer to these theories as phase-only theories. At small values of the momenta, 
the dispersion relation of the sound mode is linear in momentum, i.e., uis(k) = Cok. By using 
the ansatz but removing the breathing modes, we find the phase-only sound velocity 



(22) 



c = X^/JBoU. 

This agrees with the results of Refs. H El El El H . 

However, this result is changed when the transverse breathing mode is taken into account. 
Using the unified description of sound and breathing modes that we have presented in this 
paper, we obtain a different sound velocity, namely 
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FIG. 2: The relative shift in the sound velocity c from the phase-only result in Eq. (|22|) as a 
function of the interaction strength. 



To obtain this simple expression we ignored contributions proportional to J 2 as these are in 
general very small. 

In Fig. (J2J) we show the relative shift c/cq — 1 in the sound velocity from the phase-only 
theory. In the non-interacting limit the two results coincide, but the shift becomes more 
pronounced with increasing interaction strength. In particular, the real sound velocity is 
smaller than the one predicted by the phase-only theories. Asymptotically we obtain in the 
limit of strong interactions Hindoo (c/c — 1) = a/3/4 -1« —0.13. While this shift is not 
enormous, it is still large enough to be kept in mind when quantitative results are required. 
Intuitively, the reduction of the real sound velocity due to the transverse breathing mode is 
expected, since the transverse degree of freedom makes the condensate less stiff. Rather than 
all the energy of the density disturbance being pushed forward in the longitudinal direction, 
some of the excess energy is lost in exciting the transverse degree of freedom. 

It is interesting to observe that while a "pure" breathing mode, unaffected by the global 
phase dynamics, is possible when all the sites are breathing in phase, a "pure" sound mode 
is never possible. The sound mode relies on global phase differences and number fluctuations 
between sites, and these will inevitably couple into the transverse breathing mode. 



B. The transverse breathing mode 



At small values of momenta the transverse breathing mode in the absence of sound mode 



has a dispersion relation 
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u; = 2 + ^(b + ^)k 2 . (24) 



However, in line with our previous result for the sound mode, Eqs. ()20p and (|2*Tj) predict a 
different behaviour for the term proportional to k 2 . Namely, we have 

u B = 2 + — (b + — + ^ )k 2 , (25) 
16 V 2£ 2 J ' K ' 

where we again ignored contributions proportional to J 2 . 

The breathing-mode dispersion relations can be given in terms of an effective masses mj 
and m* B for the transverse breathing mode as u) = 2 + k 2 /2m,Q and ujb = 2 + k 2 /2m* B . 
In Fig. (j2J) we show the relative shift m B /mQ — 1 in the effective mass of the transverse 
breathing mode. The relative shift in the effective mass is always quite large and quickly 
approaches an asymptotic value 3/5 as the strength of the interaction increases. This shift, 
therefore, must be included in making quantitive predictions for upcoming experiments. 



IV. TIME EVOLUTION OF THE COUPLED SYSTEM 

Earlier in this paper we solved the equations of motion for the sound and the transverse 
breathing mode by assuming plane-wave solutions. This enabled us to obtain analytic solu- 
tions for the dispersion relations, but our theory can be used to solve also more complicated 
problems. In this section we demonstrate this by solving the coupled dynamics of the den- 
sity fluctuations and the transverse breathing mode, when the initial state of the condensate 
has a density dip. This problem is interesting since the earlier experiments on a condensate 
sound mode first created a density dip and then tracked the evolution of the condensate 
density 0, Q| . In this case calculating the time evolution of the coupled system of the sound 
mode and the transverse breathing mode is too complicated to be attacked analytically, but 
the problem can be readily tackled numerically. 

In Fig. 0] we show an example of the typical time evolution. In this figure we prepare the 
system with a Gaussian density disturbance 5 n (t = 0) = — 0.1 exp [— (n/10) 2 ] and then let 
it evolve. The density minimum splits into two parts propagating into opposite directions. 
The propagation velocity of the disturbance is in agreement with our result for the sound 
velocity in Eq. (|23|). but is different from the phase-only result in Eq. (|22|). The excess 
energy of the density dip excites the transverse breathing mode that remains well localized 
in the center of the lattice and only slowly spreads out further. This slow spreading of the 
breathing mode is hardly visible at the relatively short timescale of the figure. 

In Fig. |5] we demonstrate the "inverse" problem of an initial state that has a localized 
transverse breathing-mode disturbance and a homogeneous density distribution. The time 
evolution of the density distribution is then more complicated, but the time evolution of the 
breathing-mode amplitude is similar to that in Fig. 0] In Fig. El the fractional change of 
the condensate size was —25 % in the center of the lattice and such a deformation is quite 
large. Despite this strong deformation, the magnitude of the density disturbance remains 
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FIG. 3: The relative shift in the effective mass m* B of the transverse breathing mode from the 
uncoupled result in Eq. (|24j) as a function of the interaction strength. 



FIG. 4: The time evolution of (a) the density disturbance 5 n (t) and (b) the transverse breathing- 
mode amplitude |e n (i)| 2 . We used J = 0.1, U = 100, and the number of sites was 101. The initial 
state had a Gausssian density disturbance 5 n (t = 0) = — 0.1 exp — (n/10) 2 . In the figure dark 
color indicates the disturbance. It can be seen how the initial density dip splits into two dips 
propagating into opposite directions whereas at this timescale the breathing mode remains well 
localized around the location of the initial density disturbance. 



small, below 1 %. Therefore, under these conditions it would be difficult to image the density 
disturbance experimentally. 



V. SUMMARY AND CONCLUSIONS 

We have presented a unified theory of the sound propagation and the transverse breathing 
mode of a Bose-Einstein condensate in a one-dimensional optical lattice. Using a variational 
ansatz we calculated the dispersion relations of both modes and found out that the dispersion 
relations are quite strongly modified by the coupling between the sound and the transverse 
breathing mode. These changes are large enough that they should be included when making 
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FIG. 5: The time evolution of (a) the density disturbance S n (t) and (b) the transverse breathing- 
mode amplitude |e„(i)| 2 . We used J = 0.1, U = 100, and the number of sites was 101. The initial 
state had a Gausssian transverse breathing mode disturbance e n (t = 0) = Bq/2 exp — (ra/10) 2 . 
This choice corresponds to the fractional change of —25 % of the condensate size in the center of the 
lattice. In (a) dark color indicates a region of low density. In (b) it would seem that the breathing- 
mode amplitude is initially zero. This is a result of plotting the amplitude squared |e n (i)| 2 as 
opposed to plotting just the real part e' n (t). Typically the imaginary part of the breathing mode 
amplitude has a larger magnitude and therefore dominates in the amplitude squared plot. 



quantitative predictions for experiments. In principle, the sound mode is not only coupled 
to the transverse breathing mode, but also to modes with higher energy that have the same 
symmetry. In this paper we have ignored such higher-order effects, since the overlap with the 
transverse breathing mode is the largest and therefore the coupling to this mode dominates. 

In principle, the theories studying two weakly-coupled condensates ji^, HH H^l would 
also be influenced by the mechanism we have discussed in this paper. This makes the study 
of the phase dynamics of the weakly-coupled condensates in the presence of the transverse 
degrees of freedom an interesting topic for further research. 

When the sound mode in an ordinary elongated Bose-Einstein condensate is excited, by 
using for example a blue-detuned laser beam as in the experiments of Andrews et al. |3|,l^ ; 
the transverse profile of the condensate is modified close to the laser beam. After removing 
the laser beam, the density dip starts propagating with the sound velocity. However, we 
have seen that in principle in addition to this also a surface disturbance starts propagating 
in the condensate. This is analogous to the problem we have studied in this paper. In 
contrast with the Bose-Einstein condensate in an optical lattice, in a cigar shaped three- 
dimensional Bose-Einstein condensate the influence of the transverse degrees of freedom is 
expected to be small, and theories ignoring them p, LD, |8( are indeed in agreement with the 
experiments 0, 0, 0] • 

Our theory can also be used to study how the transverse breathing mode excites density 
modulations of the condensate in an optical lattice. In this paper we gave one numerical 
example along these lines, but when the transverse breathing-mode amplitude is very large 
or strongly modulated, we expect that nonlinearities will play an important role. Under 
such conditions instabilities might arise, revealing new possibilities for studying nonlinear 
matter wave dynamics. 
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